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1.  Introduction 

In  large  eddy  simulation  (LES)  of  turbulent  flows  the  dynamics  of  the  large 
scale  structures  are  computed,  while  the  effect  of  the  small  scale  turbulence 
is  modeled.  The  differential  equations  describing  the  space-time  evolution 
of  the  large  scale  structures  are  formally  derived  by  applying  a low-pass 
filter  with  non-uniform  filter  width  to  the  Navier-Stokes  equations.  For  an 
incompressible  flow  the  non-dimensional  equations  describing  the  evolution 
of  large  scale  structures  take  the  following  form: 
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where  the  square  bracket  denotes  the  commutation  operator  given  by 
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The  filtered  convective  term  ufuf  is  unknown  in  LES  and  is  typically  de- 
composed into  the  convective  term  ufuj  that  can  be  computed  and  the 
remainder,  called  sub-grid  scale  (SGS)  stress,  which  should  be  modelled: 
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In  order  to  derive  LES  equations  from  Eqs.  (1),  (2),  and  (4)  the  following 
three  key  assumptions  are  made: 

1.  The  differentiation  and  filtering  operations  commute. 

2.  The  shape  of  the  low-pass  filter  is  known. 

3.  The  modeled  subgrid  scale  stress  corresponds  to  the  low-pass  filter 
used  in  deriving  LES  equations. 

With  these  assumptions  we  obtain  the  classical  LES  equations 
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Ironically,  all  three  assumptions  used  in  derivation  of  Eqs.  (5)-(6)  are  very 
often  violated,  perhaps  as  a result  of  addressing  other  very  important  prob- 
lems related  to  numerical  methods.  However,  recent  advancement  in  numer- 
ical algorithms  and  hardware,  has  necessitated  a careful  and  systematic 
examination  of  this  issue.  The  main  objective  of  this  paper  is  to  present  a 
way  for  LES  to  be  consistent,  i.e.  satisfy  all  these  assumptions. 


2.  Implicit  Filtering 


Due  to  the  lack  of  a straightforward  and  robust  filtering  procedure  for 
inhomogeneous  flows,  most  large  eddy  simulations  performed  to  date  have 
not  made  use  of  explicit  filtering.  The  nearly  universal  approach  for  LES  in 
complex  geometries  is  to  argue  that  the  finite  support  of  the  computational 
mesh  together  with  the  low-pass  characteristics  of  the  discrete  differencing 
operators  effectively  act  as  a filter  (Rogallo  and  Moin,  1984;  Lund,  1997), 
i.e. 
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where  S/Sx  and  d/dx  are  respectively  finite  difference  and  the  exact  (analyt- 
ical) derivative  operators.  This  procedure  is  commonly  referred  as  implicit 
filtering  since  an  explicit  filtering  operation  never  appears  in  the  solution 
procedure.  Although  the  technique  of  implicit  filtering  has  been  used  ex- 
tensively in  the  past,  there  are  several  compelling  reasons  to  adopt  a more 
systematic  approach.  Foremost  of  these  is  the  issue  of  consistency.  While  it 
is  true  that  discrete  derivative  operators  have  a low-pass  filtering  effect,  the 
associated  filter  acts  only  in  one  spatial  direction  in  which  the  derivative  is 
taken.  This  fact  implies  that  each  term  in  the  Navier-Stokes  equations  is 
acted  on  by  a distinct  one-dimensional  filter  and  the  actual  equations  being 
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solved  are 
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where  (•)  and  (■)  ' are  the  effective  one-dimensional  filters  associated  with 
the  first  and  second  difference  operators  respectively.  Thus,  there  is  no  way 
to  derive  the  discrete  equations  through  the  application  of  a single  three- 
dimensional  filter.  Considering  this  ambiguity  in  the  definition  of  the  filter, 
it  is  nearly  impossible  to  make  detailed  comparisons  of  LES  results  with 
filtered  DNS  or  experimental  data. 

The  second  significant  limitation  of  the  implicit  filtering  approach  is 
the  inability  to  control  numerical  error.  Without  an  explicit  filter,  there 
is  no  direct  control  of  the  energy  in  the  high  frequency  portion  of  the 
spectrum.  Significant  energy  in  this  portion  of  the  spectrum  coupled  with 
the  non-linearities  in  the  Navier-Stokes  equations  can  produce  significant 
aliasing  error.  Furthermore,  all  discrete  derivative  operators  become  rather 
inaccurate  for  high  frequency  solution  components  and  this  error  interferes 
with  the  dynamics  of  the  small  scale  eddies.  This  error  can  be  particularly 
harmful  when  the  dynamic  model  is  used  (Germano  et  al.,  1991;  Ghosal 
et  al.,  1995),  since  it  relies  entirely  on  information  contained  in  the  smallest 
resolved  scales. 


d2^Xl  d2%X2  d2ST3 
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3.  Explicit  Filtering 

The  difficulties  associated  with  the  implicit  filtering  approach  can  be  al- 
leviated by  performing  an  explicit  filtering  operation  as  an  integral  part 
of  the  solution  process  (Ghosal,  1996;  Vasilyev  et  al.,  1998).  By  damping 
the  energy  in  the  high  frequency  portion  of  the  spectrum  it  is  possible  to 
reduce  or  eliminate  the  various  sources  of  numerical  error  that  dominate 
this  frequency  range.  Explicit  filtering  reduces  the  effective  resolution  of  the 
simulation,  but  allows  the  filter  size  and  shape  to  be  chosen  independently 
of  the  mesh  spacing.  In  addition  explicit  filtering  provides  the  means  for 
both  assessing  and  minimizing  the  effects  of  numerical  error  in  practical 
simulations. 

In  order  to  illustrate  the  explicit  filtering  procedure  further,  consider  a 
discrete  time  integration  applied  to  the  LES  Eqns.  (5)-(6): 

K+1  =M?  + At«J,+1  + AArL«?+1, 


(10) 
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where  A i/u”+1  and  A atlw”+1  are  incremental  velocity  fields  resulted  from 
the  time  integration  of  linear  and  non-linear  terms  respectively.  The  incre- 
mental velocity  field  A has  frequencies  beyond  the  characteristic 
frequency  that  defines  Ui  and  thus  should  be  filtered.  This  additional  filter- 
ing should  be  explicitly  performed  at  the  end  of  each  time  step.  Thus  the 
modified  time  integration  procedure  can  be  schematically  written  as 

u”+1  =u?  + A Lu?+l  + A nlu^K  (11) 

Note  that  the  frequency  content  of  each  term  on  the  right-hand  side  of 
Eq.  (11)  is  limited  to  the  bar  level.  Thus  in  advancing  from  time  level  n to 
n + 1,  the  frequency  content  of  the  solution  is  not  altered.  This  fact  implies 
that  the  additional  filtering  of  the  non-linear  terms  is  sufficient  to  achieve 
an  explicit  filtering  of  the  velocity  field  for  all  time. 

The  explicit  filtering  procedure  (11)  automatically  implies  an  alternative 
form  of  LES  equations  (Vasilyev  et  al.,  1998) 

dui  dxiiUj  _ dp  1 d2Ui 

dt  dxj  dxi  Re  dxjdxj 

where  Tij  denotes  the  alternative  SGS  stress 


T ij  — Itj  Uj  Ui  Uj  • 

A possible  drawback  of  this  new  formulation  is 
Galilean  invariant  unless  a sharp  cut-off  filter  is 
1998).  Non-Galilean  invariance  follows  from  the 
Cjd(ui  — Uj)/dxj,  where  Cj  is  the  uniform  translation  velocity  of  the  moving 
frame  of  reference.  The  error  is  seen  to  be  proportional  to  the  difference 
between  the  singly  and  doubly  filtered  velocity.  This  difference  will  be  zero 
for  a sharp  cutoff  filter,  but  will  not  vanish  in  the  general  case.  The  spectral 
content  of  the  error  is  proportional  to  G(k)(l— G(k))  where  G(k)  is  the  filter 
transfer  function  and  k is  the  wave  vector.  This  fact  implies  that  the  error 
is  only  generated  in  the  wavenumber  band  where  G(k)  differs  significantly 
from  0 or  1.  Thus  it  is  possible  to  minimize  the  error  by  constructing  the 
explicit  filter  to  be  as  close  as  possible  to  a sharp  cut-off. 

It  turns  out  that  non-Galilean  invariance  of  Eq.  (12)  is  not  an  issue  if 
one  adapts  the  correct  physical  interpretation  of  inhomogeneous  turbulent 
complex  geometry  flows.  We  recall  that  if  the  filter  width  is  non-uniform 
in  one  or  more  directions,  it  automatically  implies  that  the  LES  equations 
should  not  be  Galilean  invariant  in  these  directions  (Ghosal,  1999).  To  en- 
sure Galilean  invariance  in  homogeneous  directions  of  the  flow  one  may  use 
sharp  cut-off  filters  in  these  directions. 
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that  the  Eq.  (12)  is  not 
utilized  (Vasilyev  et  al, 
appearance  of  the  term 
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Figure  1.  Modified  wave- number  diagram  for  the  second  (left)  and  fourth  (right)  order 
schemes.  The  vertical  solid  line  represents  the  effective  resolution  of  the  explicitly  filtered 
simulations,  i.e.  fixed  cut-off  frequency  k.  The  dashed  line  is  the  exact  distribution  that 
is  achieved  with  a pseudo  spectral  method. 


It  is  worthwhile  to  note  that  the  use  of  explicit  filtering  provides  the 
means  for  reducing  the  various  sources  of  numerical  error  that  become 
most  severe  for  length  scales  on  the  order  of  the  mesh  size.  By  damping  the 
high  frequency  portion  of  the  solution,  it  is  possible  to  control  the  adverse 
effects  of  numerical  error.  In  particular,  if  the  filter  width  is  held  fixed  as 
the  mesh  is  refined,  the  velocity  field  will  converge  to  the  true  solution  of 
LES  equations  (5)  and  (12).  This  should  be  contrasted  to  the  conventional 
approach  where  the  mesh  is  refined  without  the  use  of  an  explicit  filter. 
In  the  latter  case,  additional  length  scales  are  added  each  time  the  mesh 
is  refined,  and  thus  the  process  converges  to  a direct  numerical  simulation 
(DNS)  rather  than  to  an  LES. 

It  is  very  illustrative  to  consider  the  effect  of  explicit  filtering  on  the 
accuracy  of  the  solution  and  the  computational  cost  associated  with  it 
in  terms  of  modified  wavenumber  analysis  (Lund  and  Kaltenbach,  1995). 
For  simplicity  let  us  consider  the  case  when  we  filter  the  periodic  signal 
given  on  a uniform  mesh  with  grid  size  Ag.  Then  explicit  filtering  can  be 
viewed  as  refining  the  mesh  but  keeping  the  same  frequency  content  ( k < 
k)  determined  by  the  primary  filter  width  A.  The  modified  wavenumber 
diagrams  for  these  simulations  are  shown  in  Fig.l.  The  solid  vertical  line 
denotes  the  fixed  effective  resolution,  while  the  curves  to  the  left  of  this 
line  show  the  modified  wavenumber  distributions  for  the  various  values 
of  filter  width  to  grid  ratios  A/Ag.  When  no  filter  is  applied  (the  case 
A/Ag  = 1)  considerable  truncation  error  is  evident  for  the  upper  half 
of  the  wavenumber  range.  As  the  ratio  A/Ag  is  increased,  the  situation 
improves  and  the  modified  wavenumber  slowly  approaches  the  analytical 
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distribution  for  k < k.  Note  that  with  the  increase  of  the  ratio  A/Ag  the 
computational  cost  associated  with  the  implementation  of  explicit  filtering 
also  increases.  Thus,  in  order  for  explicit  filtering  to  be  cost  effective,  the 
higher  order  schemes  need  to  be  used,  since  they  provide  the  comparable 
modified  wavenumber  dependence  in  the  range  k < k with  considerably 
smaller  filter  width  to  grid  ratios  then  the  lower  order  schemes. 

In  summary,  the  use  of  explicit  filtering  results  in  a LES  that  satisfies 
assumptions  2 and  3 that  are  used  in  the  derivation  of  the  LES  equations. 
However,  in  order  to  realize  the  benefits  of  explicit  filtering  and  satisfy 
assumption  1,  it  is  necessary  to  develop  discrete  filtering  operators  that 
commute  with  numerical  differentiation.  This  will  be  discussed  next. 


4.  Construction  of  Commutative  Filters 


In  order  to  realize  the  benefits  of  explicit  filtering,  it  is  necessary  to  develop 
robust  and  straightforward  discrete  filtering  operators  that  commute  with 
numerical  differentiation.  Recently  a new  class  of  commutative  filters  for 
both  structured  (Vasilyev  et  al.,  1998)  and  unstructured  (Marsden  et  al. , 
2000)  grids  has  been  developed.  With  these  filters  the  differentiation  and 
filtering  operations  commute  to  an  a priori  specified  order  of  filter  width. 

We  begin  by  discussing  the  continuous  filtering  operation  and  then  ex- 
tend it  to  discrete  filtering.  Let  us  consider  a three  dimensional  field  <£(x) 
(x  = (xi,  X2,  xa)T)  defined  in  a three-dimensional  domain  D.  The  most  gen- 
eral continuous  filtering  operation  in  physical  space  can  be  written  as 

<?(x) = ^ G (x>  y)  ^(y)rf3y>  (14) 


where  G (x,  y)  is  the  location  dependent  three-dimensional  filter  function. 
The  filtering  operation  (14)  can  be  substantially  simplified  if  there  exists 
a function  £ = f(x)  and  its  inverse  x = F(£)  that  maps  the  physical 
domain  D into  a rectangular  domain  = [a\,f}\]  x [as,  A>]  x [03,^3].  In 
this  case  the  filtering  operation  can  be  defined  the  following  way:  given 
an  arbitrary  function  <f>(x)  we  obtain  the  new  function  rp(^)  = </>(F(£)) 
defined  in  the  domain  fi,  the  function  ■!/>(£)  is  then  filtered  using  a sequence 
of  three  one-dimensional  filters.  It  can  be  shown  (Vasilyev  et  al.,  1998)  that 
the  commutation  error  of  filtering  and  differentiation  operations  in  three 
spatial  dimensions  is  given  by 


d<j>  _ d<f>  d<f> 

dxk\  ~ dxk  dxk 


0(A”,  A2,  A3), 


(15) 


provided  that  the  filter  moments  M*(£)  of  filters  Gi  (i  = 1, . . . , 3)  exist  and 
they  satisfy  the  following  properties: 
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Figure  2.  (a)  Fourier  transform  G(k)  of  the  symmetric  discrete  filters  with  three  van- 

ishing moments  and  different  linear  constraints,  (b)  Real,  imaginary,  and  absolute  value 
of  G(k)  of  the  asymmetric  discrete  filter  with  three  vanishing  moments. 

M?(6)  = 1 for  & e [a A],  (16) 

M*(£i)  = 0 for  k = 1,. . . ,n  — 1 and  & € [a*,  A].  (17) 

Also  it  can  be  shown  that  the  commutation  error  in  the  homogeneous  di- 
rections is  exactly  zero. 

The  discrete  filtering  operation  is  defined  as 

h 

0,  = S w/ <£?+/’  (18) 

l=—Kj 

where  wj  are  the  filter  coefficients  chosen  to  satisfy  the  zero-moment  re- 
quirements (16),  (17)  and  possibly  additional  constraints.  Examples  of  dis- 
crete symmetric  and  asymmetric  filters  are  shown  in  Fig.  2. 

The  discrete  commutative  filter  construction  theory  can  be  easily  ex- 
tended to  the  case  when  mapped  space  does  not  exist,  e.g.  when  unstruc- 
tured meshes  are  used.  In  this  case  there  is  no  choice,  but  to  construct 
filters  in  the  physical  domain.  In  this  case  it  can  be  shown  (Marsden  et  al, 
2000)  that  for  a smoothly  varying  filter  width,  the  local  commutation  error 
in  three  dimensions  is  given  by 

~ =0( Al1(x)AJ2(x)A^(x))  , i+j  + k = n (19) 

provided  that  filters  moments  MIjfc(x)  = — fq,  ( T] , x)  d3rj  exist  and 

satisfy  the  following  constrains: 

= o j<i7i  + k< 


(20) 
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Figure  3.  Example  of  filter  constructed  with  triangles  on  an  unstructured  grid  and  the 
corresponding  transfer  function. 


Greater  predictability  and  ease  of  implementation  can  be  gained  by  us- 
ing interpolation  based  filters  that  are  constructed  by  evaluating  either  in- 
terpolating polynomials  (Marsden  et  al.,  2000)  or  least  square  interpolants 
(Haselbacher,  2001)  at  the  desired  locations.  The  nth  polynomial  construc- 
tion automatically  guarantees  an  nth  order  commutation  error.  To  control 
the  shape  and  other  properties  of  the  discrete  filters,  filters  can  be  con- 
structed as  a linear  combination  of  many  polynomial  based  filters,  while 
preserving  the  commutation  properties  of  the  filter. 

In  general,  with  an  nth  order  numerical  scheme,  the  filtering  operation 
must  commute  to  order  n.  Most  of  the  unstructured  mesh  codes  are  sec- 
ond order  accurate,  thus  a discrete  filter  based  on  first  order  polynomial 
interpolation  is  sufficient  to  achieve  second  order  commutation  error.  A 
two-dimensional  discrete  filter  based  on  first  order  polynomial  interpola- 
tion can  be  constructed  using  a triangle  of  three  points  surrounding  the 
point  (xo,yo)  where  we  want  the  filtered  value.  A triangle  is  chosen  be- 
cause in  two  dimensions  three  points  are  needed  for  exact  reconstruction  of 
a first  order  polynomial.  Weights  are  calculated  by  fitting  a polynomial  to 
the  vertices  of  the  triangle,  and  are  then  used  to  find  a weighted  average  at 
the  central  point  (xq,  yo).  With  this  method,  the  same  number  of  points  are 
used  in  the  filter  for  any  point  on  the  mesh.  Figure  3 shows  an  example  of 
a 2-D  filter  constructed  from  three  triangles.  To  achieve  flexibility  in  filter 
width,  each  triangle  as  well  as  the  central  point  is  assigned  a weight  Pi, 
which  applies  equally  to  all  vertices  of  the  triangle.  The  value  of  Pi  can  be 
varied  from  0 to  1 as  long  as  the  sum  total  is  1.  The  optimum  value  of  p 
for  the  central  point  is  1/2  because  this  results  in  a transfer  function  with 
a well  defined  peak  and  low  pass  filter  shape. 
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5.  Spectrum  of  Commutation  Error 

Leading  order  commutation  error  analysis  described  in  the  previous  section 
can  be  viewed  as  a practical  tool  for  constructing  discrete  filters  that  com- 
mute with  finite  difference  operators  to  an  a priori  specified  order  of  filter 
width.  However,  the  leading  order  error  analysis  by  itself  is  not  enough  to 
guarantee  that  the  commutation  error  is  negligible  compared  to  the  subgrid 
scale  stress,  since  it  does  not  use  any  information  about  spectral  content 
of  the  analyzed  signal  (Ghosal,  1996;  Ghosal,  1999).  Due  to  the  presence 
of  significant  energy  in  the  high  frequency  portion  of  the  LES  spectrum, 
the  commutation  error  could  be  considerable  and  in  some  cases  even  com- 
parable with  the  subgrid  scale  stresses.  The  existence  of  inhomogeneous 
directions  in  the  flow  precludes  the  use  of  the  standard  Fourier  analysis. 
In  this  section  we  present  a recently  developed  mathematical  tool  that  can 
be  used  to  analyze  the  local  spectrum  of  the  commutation  error  and  its 
dependence  on  the  filter  shape  and  the  non-uniformity  of  the  filter  width. 

We  begin  by  introducing  the  local  spectrum  analysis  in  one-dimensional 
space  and  then  extend  it  to  three  spatial  dimensions.  Let  us  consider  a one 
dimensional  filter  of  constant  shape  but  variable  width.  The  continuous 
filtering  operation  can  be  written  as 

?w=a hflG(m)'KvW  ,21) 

Substituting  the  Fourier  integral 

r $ (Jfc)  eiKXdk  (22) 

27T  J- oo 

into  Eq.  (21),  changing  the  order  of  integration,  and  integrating  the  result- 
ing equation  with  respect  to  y we  obtain  the  following  equation: 

<f>{x)  = — G[A(x)k)4>{k)elkxdk,  (23) 

27 r J — oo  v ' 

where  G ( k ) is  the  Fourier  transform  of  the  filter  function  given  by 

/OO 

GWe-^dt.  (24) 

-OO 

Now  comparing  Eqs.  (22)  and  (23)  we  can  see  that  the  structure  of  the 
equations  is  the  same.  Thus,  the  local  Fourier  transform  can  be  defined  as 


(f)(k,  x)  — G (A  (x)  k ) 4>(k). 


(25) 
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G(k)  K(k) 


Figure  4-  Transfer  functions  G(k ) and  K(n)  for  a variety  of  filters.  Mn  denotes  the 
discrete  filter  with  n vanishing  moments. 


Now  the  meaning  of  the  local  Fourier  transform  is  clear,  it  refers  to  the 
location  in  space  where  the  transform  is  taken  and  reflects  the  fact  that 
the  filter  width  is  a function  of  location.  Performing  analogous  analysis  for 
commutation  error  it  can  be  shown  that 


dcf)' 

dx 


A'(x)  1 
A(x)  27r 


/o° 

K (k)  <fi(k)elkxdk, 

-OO 


(26) 


where  k = A {x)  k and  the  transfer  function  K (k)  is  defined  by 


K(k)  = -k^^.  (27) 

diKi 

Now  the  effect  of  filtering  is  clearly  seen:  filter  width  stretching  affects  the 
amplitude  of  the  commutation  error,  while  the  filter  shape  and  width  affect 
the  spectrum.  The  effect  of  the  filter  shape  on  the  spectrum  of  commutation 
error  is  demonstrated  in  Fig.  4,  where  transfer  functions  G(k)  and  K (k) 
are  shown  for  a variety  of  filters.  It  is  important  to  note  that  the  closer  the 
filter  to  the  sharp  cut-off,  the  more  localized  the  spectrum  of  commutation 
error.  Also  note  that  the  commutation  error  is  exactly  zero,  when  the  width 
of  the  filter  is  constant  throughout  the  domain. 

The  one-dimensional  analysis  can  be  easily  extended  to  three  spatial 
dimensions.  In  particular,  it  can  be  shown  that  for  three  dimensional  filters 
of  constant  shape  and  variable  width  the  commutation  error  is  given  by 


d4> ' 

dxi . 


1 dA;(x)  1 
A;(x)  dxi  (2tt)3  7-oo 


Ki  (k)  <^>(k)eikxd3k, 


(28) 


where  A j(x)  is  the  filter  width  in  the  Xi  direction,  k = (k\,  ^3)T  and 

K = (Ai(x)fci, A2(x)&2,  A3(x)/c3)t  are  three  dimensional  wave  vectors, 
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Figure  5.  One-dimensional  spectrum  of  the  exact  SGS  force,  drn/dxi , and  local 
one-dimensional  spectrum  of  the  commutation  error  for  forced  isotropic  turbulence  at 
Re\  = 168  for  (a)  Gaussian  and  (b)  My  filters. 


4>{ k)  is  the  Fourier  transform  of  function  </;(x),  G(k)  is  the  Fourier  transform 
of  the  filter  function,  and  the  transfer  function  Ki  (k)  is  defined  by 


Ki  (k)  = -Kj 


0G(k) 

Dki 


(29) 


To  demonstrate  the  effect  of  commutation  error  in  LES,  we  consider 
the  results  of  2563  direct  numerical  simulation  of  forced  homogeneous  tur- 
bulence at  Re\  = 168  (Jimenez  and  Wray,  1993).  Figure  5 shows  the 
one-dimensional  spectrum  of  the  exact  SGS  force,  drn/dxi,  and  the  lo- 
cal one-dimensional  spectrum  of  commutation  error  for  the  Gaussian  and 
the  discrete  M7  filter.  These  two  filters  are  chosen  because  the  Gaussian 
filter  is  smooth,  while  the  M-j  filter  is  close  to  the  sharp  cut-off  filter.  The 
exact  SGS  force  is  obtained  by  filtering  the  DNS  data.  Note  that  the  lo- 
cal spectrum  of  the  commutation  error  shown  in  Fig.  5 does  not  take  into 
account  the  grid  stretching  factor.  The  results  confirm  strong  dependence 
of  the  spectrum  of  commutation  error  on  the  filter  shape:  the  spectrum  is 
global  for  smooth  filters  like  Gaussian,  while  for  filters,  that  are  close  to 
the  sharp  cut-off,  the  spectrum  is  localized. 


6.  Conclusions 

In  this  paper  we  presented  a systematic  approach  for  large  eddy  simulation 
to  be  consistent,  i.e.  to  satisfy  all  the  assumptions  used  in  the  derivation 
of  the  LES  equations.  We  demonstrated  that  a solution,  consistent  with 
the  true  filtered  Navier-Stokes  equations,  can  be  achieved  by  means  of  ex- 
plicit filtering.  We  have  shown  that  explicit  filtering  not  only  provides  the 
means  for  both  assessing  and  minimizing  the  effects  of  numerical  error  in 
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practical  simulations,  but  also  allows  the  filter  shape  and  size  to  be  chosen 
independently  of  the  computational  mesh.  Having  the  ability  to  control  the 
filter  shape  opens  whole  new  horizons  in  LES  modelling,  since  it  makes 
it  feasible  to  look  at  SGS  modelling  and  filtering  as  one  inseparable  issue. 
The  success  of  explicit  filtering  strongly  depends  on  the  ability  to  construct 
discrete  filters  that  commute  with  differentiation.  We  have  described  the 
general  theory  of  construction  of  such  filters  and  developed  a tool  for  as- 
sessing the  local  spectrum  of  the  commutation  error. 

Although  tremendous  progress  has  recently  been  made  in  the  develop- 
ment of  the  explicit  filtering  approach,  a lot  of  fundamental  research  needs 
to  be  done  before  researchers  and  engineers  can  take  full  advantage  of  this 
new  mathematical  tool.  Among  the  issues  which  definitely  need  to  be  ad- 
dressed in  the  near  future  are  the  development  of  the  SGS  models  that 
correspond  to  specific  filters  and  their  subsequent  verification  with  experi- 
mental and  DNS  results. 
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